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PREFACE 


Knowledge of thermal conditions in the ground is necessary for the planning 
and performing of various kinds of building work. In house-building activities, 
the depth of the permanently unfrozen ground (the course of the 0°C isotherm) 
has to be known; this applies also to road, aerodrome and other construction 
work. Various kinds of pipelines, such as those for water mains, sewerage and 
drainage, have to be laid out at a depth where the ground will not freeze. Cal- 
culation of the heat losses in remote heating ducts requires knowledge of the ther- 
mal conditions in the ground, Heating systems based on the heat pump principle 
mostly employ the heat stored in the ground for their heat source, 

The present investigation was instigated by a study of the floor heating of build- 
ings erected upon a concrete slab poured directly on the ground, which the au- 
thor has carried out in the State Institute for Technical Research in Finland. 
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1, FACTORS AFFECTING THE TEMPERATURE OF THE UPPER GROUND LAYERS 


The field of temperature distribution in the ground is determined by the phys- 
ical quality of the ground and by the entering and escaping quantities of heat. 
The ground obtains its entire heat from the sun by direct solar radiation (S,), dif 
fuse, scattered radiation from the atmosphere (S},) or counterradiation from the 
atmosphere (Sy). Owing to the geothermal gradient, about 0,03°C/m in magni- 
tude, a quantity of heat amounting only to about 0,06 kcal/m%h flows from the 
interior of the earth to its surface [1]. This heat quantity is negligible with re- 
spect to its effect upon the temperature of the upper ground layers. 

On the other hand the ground immediately loses part of the received radiant 
energy in the form of radiation reflected by the soil surface (S;). According to 
Angstrom’ s [2] and Lunelund’s [3] investigations, the proportion of reflected ra- 
diation has the following magnitude, depending on the quality and covering (veg- 
etation; snow cover) of the ground: 


Quality of soil surface Reflected quantity of heat, 
per cent, (Albedo) 
Grey sandy plain 12 26 
Black earth 12“ 14 
Granite rock 1218 
Cereal fields, on an average 16 
Grass plain 25» 26 
Forest 10 18 
Pure snow cover 81 85 
Old snow cover 42-70 


The reflected quantities of heat vary greatly according to the quality of the 
ground, 
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Part of the radiant energy retained by the ground is continuously given off as 
radiation into space; the soil surface acts approximately in the manner of a black 
body giving off long-wave thermal radiation (S,). 

The difference of the heat quantities received and given off by the ground, 
the radiation balance (S,), determines the temperatures of the upper ground layers 
(010 m) and the magnitude of the heat quantities stored in the ground. 


St = (Sg + Sp + Sy) = (Sp + Sy) 


1 
= (Sy + S, (S,- 


The annual quantities of radiant energy in Finland are, according to Kerdnen 


[4]: 


S, = 446 Mcal/m? S, + S, = S, = 640 Mcal/m” 
331" 
211" 


During the warm season, the ground receives more heat in the form of radi- 
ation than it gives off during the whole 24-hour day, i.e., the heat balance (S;) 
of the ground is positive and the temperature of the surface layers increases. 

According to Keranen [4], the monthly thermal balance values in South-Finland 
are: 


Sy I | V | Vi | IX | X | XL | XM 


Mcal/m? | «29 | -23 | -2 | 26 | 61 | 95 | 82 | 47 | 16 }-12|-23]-27 | 211° 


The thermal balance is positive during the time from April to September. 

The quantity of radiant energy that remains available on the ground surface 
(Sy) is consumed to heat the upper ground layers (Q,,), to heat the adjacent air 
(Qj) and to evaporate water and to melt ice Q, [5], [6], [7]: 


S, = Qn + Q + Q (2) 


Heating of the ground itself requires only a very small fraction of the heat 
quantities available on the soil surface. According to Franssila’s investigations [8], 
only about 5 % of the radiation balance on a warm summer day remain in the 
ground. 

A remarkable part of this radiant energy is transferred to the air by the ef- 
fects of turbulence and by conduction. Great quantities of heat are consumed to 
evaporate water, particularly in windy weather [9]. On the other hand the ground 
receives heat from warm air masses and from water condensing on its surface, re- 
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leasing its latent heat, Rain will either heat or cool the ground, depending on the 
temperature of the rain water. In Finland, rain usually increases the ground tem- 
perature. 

According to the investigations carried out by Singer [10] in Munich, the tem- 
peratures in the upper ground layers deviate upward from the mean the more, the 
greater the amount of rainfall during the preceding summer. This is mainly due 
to increase in thermal conductivity of the soil with increasing moisture. According 
to measurements by Homen [11] at Mustiala, the temperature of the surface lay- 
ers of the ground is 0,5 3,0°C higher in rainy summers than in dry summers. 

The annual range of variation of the ground temperatures and the magnitude of 
the heat quantity stored in the ground depend further on the soil quality, its mois- 
ture conditions, density, specific heat, etc, The heat transfer from the surface 
layers to greater depths is controlled by the thermal conductivity of the soil and 
by the temperature gradient in the ground, The heat stored in the ground is at 
its maximum in the beginning of September, when the heating period of buildings 
begins, and it is lowest in the beginning of April [12]. In the intervening period 
the ground gives off the heat stored during the summer. Owing to its high thermal 
insulating capacity, snow cover delays the cooling of the ground. During winters 
abounding in snow the ground remains unfrozen up to its surface, owing to the heat 
stored in the ground. 
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2 THERMAL CONDUCTIVITY AND TEMPERATURE CONDUCTIVITY OF THE 
SOIL, AND THEIR DETERMINATION 


For the mathematical calculation of the temperature distribution in the ground 


it is necessary to know a physical constant of the soil substance, its temperature 
conductivity a: 


A 
er (3) 


This quantity is dependent on three factors: the thermal conductivity A, the vol- 
ume weight@ , and the specific heat c. 

According to Kersten’s investigations, the thermal conductivity of the soil is 
dependent on its moisture @ (in percent. of the dry weight) and its dry volume 
weight @,, by the formula, 


- for coarse soil types: 


A= (0,087 log p+ 0,05) + 10°*8% + @o 


- for fine soil types: (4) 
A= (0, 115 log 0, 025) 19%» 624 . 


The constants in the formulas have been determined by the said author on the 
basis of .a material comprising several thousand A determinations under varying 
conditions of moisture, volume weight and temperature. The formulas were cal- 
culated for a mean temperature of +4,5°C and the maximum deviation of the 
observed values from those calculated by the formulas were within 25 %, when 
1 %< P< 25 % and 1400 kg/m*< < 2200 kg/m® [13]. 

As a tule, the dry volume weights of the different soil types vary within nar- 
row limits, according to the investigations of Watzinger, Kersten and Saare {13}. 
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Most soil types in the dry state have, according to Pfaundler [14], a specific 
heat, which is very accurately equal to 0,2 kcal/kg °C, or one fifth of the spe- 
cific heat of water, if the soil is free of humus substances. In highly humous 
soils, the specific heat may be as high as 0,4 kcal/kg %G, 

The heat content of moist soil can be computed from the specific heats of the 
dry soil and of water by the following formula: 


Cy = (@- wy Ce +we Cy (5) 


Heat content of the moist soil [kcal/m® 
Volume weight of the moist soil [kg/m] 

w Water content of the soil [kg/m*] 

C, Specific heat of the dry soil [kcal/kg mt 

Cy Specific heat of water [kcal/kg °C] 


In Table 1, values of the dry volume weight, moisture contént, thermal con- 
ductivity, heat capacity and temperature conductivity of the soil types occurring 
most commonly on building sites have been compiled, arranged by the tempera- 
ture conductivity of the soil. The values are results from measurements carried 
out by Kersten in Alaska, E, Saare and C.G, Wenner in Sweden and Watzinger in 
Norway. In their measurements, Saare and Wenner have employed the non-station- 
ary probe method, the soil being in undisturbed natural state. Watzinger has meas- 
ured values of the thermal conductivity by a stationary method, the soil under in- 
vestigation being placed between a heated and a refrigerated plate and subjected 
to load. Kersten’s method is otherwise the same but in his apparatus a cylinder 
replaces the heated plate and another, annular cylinder the cooled plate [13]. 
The values obtained by the non-stationary method in the natural state are as a 
rule about 10 % in excess of those obtained by the stationary method. A family 
of graphs plotted on the basis of these values, showing the dependence of the tem- 
perature conductivity of the soil types in question on their moisture, is presenicd 
in Fig. 1. 

In the nature, temperature conductivity is not a constant even if one and the 
same soil type is considered but it varies according to the moisture content and, 
in some degree, also to the temperature of the soil. At low moisture contents, a 
relatively small change in water content as much as doubles the temperature con- 
ductivity. This has also been shown by Keen and Patton [9] in their investiga- 
tions. 

The moisture in the ground varies considerably with the seasons and according 
to the abundance of rainfall, with consequent variations in temperature conductiv- 
ity is highest in the spring and autumn, when the ground is wet. In the summer, 
when the ground is usually dry, it has a low temperature conductivity. 
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Fig. 1 Thermal conductivity of different soil types in dependence of mois- 
ture content (% of dry weight), 


Key to the curves: 


Gravel (dry volume weight 1760 t 90 kg/m’), Kersten, Alas 


" (1920 40 kg/m*), 4 

" (2055 75 kg/m’), 
Sandy gravel (1800 kg/m®), Watzinger, Norway 
Natural esker gravel (1550 kg/m’), 


Gravelly medium coarse sand (1716 * 16 kg/m®), Kersten, Alas 
" "(1800 20 kg/m, =" 

Fine sand (1792 * 8 kg/m’), . ° 

" (1680 30 kg/m), 

" (1580 kg/m), Watzinger, Norway 
Very fine sand (1763 * 7 kg/m’), Kersten, Alaska 
Fine sandy silt (1462 kg/m’), 
Very heavy clay (1498 * 8 kg/m), 


" "(1337 + 50 kg/m, " 


12 

- 

my 

0006 

— 

0003 — 

' 

— 

— - 

— 

qoor 

ag 

" 

. 

‘ 


13 


Sandy and fine sandy till (2116 


36 kg/m), Kersten, Minnesota 
10 kg/m’), . 


1+ t+ 


In addition to moisture, also the annual variation of the temperature in the 
soil affects its temperature conductivity. If one excludes the surface layer of 0,5 m 
depth, the soil temperature in South-Finland may vary only within the limits of 
-5 +16°C (Figs. 5, 6 and 7), At the upper limit of this temperature range the 
thermal conductivity of the soil is at its most 5 % in excess of that at the lower 
limit. The temperature conductivity increases with increasing depth because the 
density of the soil increases. At greater depth, where the moisture is less variable 
than in the upper layers and the ground does not freeze, more reliable values are 
obtained on determining the temperature conductivity. As a result of the numerous 
factors affecting the thermal conductivity and temperature conductivity of the soil, 
calculations performed with a constant value for the temperature conductivity are 
only approximate. 

When the ground freezes, its thermal conductivity increases, as the thermal 
conductivity of ice, A = 1,94 kcal/m h °C, is considerably higher than that of 
water, A = 0,504 kcal/m h . However, coarse soil becomes more porous when 
it freezes, and this in turn reduces the thermal conductivity. For this reason the 
thermal conductivity of moist soil (below capillary saturation) changes but slightly 
when the ground freezes. In the state of capillary saturation the thermal conduc- 
tivity of the soil (not clay) increases by 10 ~ 15 % and sometimes even by 20 ~ 30%, 
when the ground freezes completely [15]. 

The heat content of ice, 0,45 kcal/dm®, is less than half the quantity of heat 
held by the same volume of water. The heat capacity will therefore decrease when 
the soil freezes; since at the same time the thermal conductivity increases, the 
temperature conductivity will also increase. The increase in temperature conduce 
tivity caused by freezing does not attain remarkable magnitude until in very wet 
soil containing more than 20 % water [16]. 
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3. THE TEMPERATURE FIELD IN FREE SOIL 


a) Calculation of the temperature field 


The temperature of free soil varies periodically in accordance with the seasons; 
the daily temperature variation is superimposed upon the annual variation. In this 
connection we are only concerned with the annual temperature wave, 

Within a small area the soil surface may be considered to be a plane, The 
isotherms in the ground are then planes as well, and the heat flow is perpendic- 
ular to them. As there is no heat flow in the lateral direction, the function rep- 
resenting the temperature distribution in the ground, T(z,t), will satisfy a one-di- 
mensional differential equation 


6 
dt dz2 (6) 


The annual temperature variation at soil surface, T(0,t) (boundary condition), 
can be represented by means of a Fourier series consisting of sine and cosine 
terms. 


co 
T@,t) = > An cosna@t + By sin nwt, (1) 


the 
: 
where 
22% 
@ = — and 
vo 
= 365 days 
Agee 
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The temperature at different depths in the ground is found by solving the dif- 
ferential equation (6) under the boundary condition (7): 


n@ 


sin (8) 


2a Zz 
T(z,t) = cos nat -- )+Bie 
2 n VYianw n 


The coefficients A, and B, are determined by the boundary condition T(0,t) = 
T(0,t) pax When t = 0. 

For the coefficient A, in equation (8) let us write A, = T, (annual mean tem- 
perature), Equation (8) becomes then 


T(z,t) = To + 
1 


e [+s cos nwt + By sinnw(t | (9) 


The annual variation is thus obtained as the sum of a number of temperature waves. 
The amplitude of each wave becomes less with increasing depth and it is subject to 
phase lag at the same time. Each wave has its own damping factor, wavelength 
and velocity of propagation in the ground. 


Wavelength: L, = ad (10) 


Velocity of propagation: v = V2anw (11) 


and 


The daily wave penetrates into the ground at higher velocity than the annual 
wave, which has a lower frequency. On the other hand, the daily temperature 
wave dies away very rapidly, extending to about 60 cm depth; it can thus be 
neglected when temperatures in ground layers at greater depth are being deter- 
mined, 

In order to calculate the annual variation of the temperatures in the ground it 
is necessary to know the temperature conductivity a of the soil type in question 
and the temperature function at soil surface, T(0,t). 
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b) The temperature function at soil surface 


No adequate data are available conceming the annual variation of the soil sur- 
face temperature but the air temperature at 2 m height from the soil surface is 
accurately known on the basis of meteorological observations in about 150 locali- 
ties in Finland [17]. The daily means of the air temperature have therefore been 
used in the subsequent calculations instead of the daily mean soil surface temper- 
atures. 

Lite investigation has been carried out in elucidation of the question how closely 
associated the variation of air temperature is with the soil surface temperatures. 


It is also a fact that the soil temperature at its very surface is exceedingly dif- 
ficult to measure. 


W. Kreutz and M. Rohweder [18] have studied this question by measuring the 
soil temperatures at 5, 10, 15 and 20 cm depth and calculating on their basis the 
soil surface temperatures, According to their experiments, the temperature of 
sandy soil follows the air temperatures measured at 5cm height within about s%. 
In Finland, measurements of temperatures at varying depth in the soil have been 
carried out by the observatory at Jokioinen in 1958+ 1959, Fig. 2 shows the an- 
nual variation of the soil temperature by daily means at 1 cm depth and the cor- 
responding mean air temperatures, as measured in a meteorological hut at 2 m 
height above soil surface, As a rule, the temperature measured at 1 cm depth in 
the soil can be considered as representing the soil surface temperature within a 
few tenths of 1°C. The temperature graphs reveal that the daily mean soil sur- 
face and air temperatures parallel each other very well. However, the soil tem- 
perature is usually on an average 1°C higher than the air temperature during the 
snow-free period. In the spring and early summer, when temperature is ascending, 
the daily mean of the soil surface temperature may rise 3 to 4°C above the mean 
air temperature, while in the autumn, when the temperature is on the downgrade, 
the soil surface temperature is at times 2 to 3°C lower than the daily mean of 
the air temperature, In the winter, when the snow forms a heat-insulating layer 
upon the ground, the soil temperature remains close to the freezing point, within 
1°C. If the soil surface is kept free of snow, the daily mean of the soil surface 
temperature follows that of the air temperature also in wintertime. 

Fig. 3 shows the annual variation of the daily mean air temperatures in Helsin- 
ki, as means of the ten-year period 1945 1954 (Curve No, 1}. Immediate Fourier 
analysis of this curve is rendered inaccurate by the sudden jumps occurring in the 
curve. It was therefore smoothed by computing ten-day means and constructing 
from them the final temperature graph (Curve No. 2), using the three-point ad- 
justing method. The mathematical equation of the curve was determined by Fourier 
analysis, determining the origin so that at the time 17. VII, 0,00 hours t = 0 and 
T(0,t) = Tray + 5,48°C, because the annual mean temperature of the ten-year 
period in question was 5,48°C, In the coordinate system defined in this manner, 
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the variation of the mean air temperature in terms of time is thus represented by 
the function 
T(0,t) = 5,48 + 11,68 cos 0,000717 t + 0,90 cos 0,001434 t + 
0,06 cos 0,00216 t» + 1,29 sin 0,000717 t - (12) 
1,04 sin 0,001434 t + 0,85 sin 0,00215 t» [ °C. ] 


c) The temperatures of free soil at various depths 


When soil temperatures are being calculated for technical purposes, the func- 
tion representing the annual variation of the soil surface temperature can be re- 
placed, with sufficient accuracy, with the temperature graph plotted from the 
daily mean temperatures obtained from meteorological observations. The calcula- 
tions will then always remain on the safe side, e.g. when the limit of ground frost 
and the heat losses are calculated, because at the most the soil temperatures ob- 
tained in this way will be 1°C too low. 

The temperatures at various depths under bare ground have been calculated from 
the temperature curve T(0,t) measured for the soil surface, by equation (13) 


T(z,t) = T, +e (13) 


The phase shift and damping coefficient has been chosen to conform with the 
first cosine wave; this wave has an amplitude about ten times that of the other 
waves and the latter die away at a much higher rate, 

The temperature distribution field has been determined for three different val- 
ues of a (Pigs. 4, 5 and 6). 


Soil type group I (Table 1) a = 0,0022 m2/h 
Soil type group Il a = 0,0032 m2/h 
Soil type group III a = 0,0047 m2/h 


When the coefficients A, and B, in equation (9) are determined by harmonic 
analysis, the temperature function T(z,t) [°C] becomes 
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T(z,t) = 6,48 + 11,68 


cos 0,000717 ( t- 
+ 1,29 0,000717 ) (14) 
1,04 ‘ ain 0,001494 (t= 


where the time t has been expressed in hours and t = 0 corresponds to July 17th, 


00,00 hours, The constants occurring in this formula have been tabulated in Ta- 
ble 2, 


Table 2 
l 
h/m 
0, 0022 0,404 563 
0, 0032 0,335 467 
0, 0047 0,275 384 


On the soil surface (z = 0) the deviation of the temperature values computed 
by this formula from the annual, adjusted mean value curve is no more than 1°C 
at any point (Figs. 4, 5 and 6). 

The amplitude of the annual temperature curve becomes less with increasing 
depth, and progressive phase shift of the wave occurs at the same time. The ve- 
locity of propagation v, wavelength L, and the depth z at which the temperature 
variation range is 1°C are seen in Table 3 for the above-mentioned soil types. 
For instance, in soil type group I (Fig. 4) the amplitude, at soil surface 12,2°C, 
is 8,2°C at 1 m depth and only 1°C at 6 m depth, The temperature extremes 
occur about 30 days later at 1 m depth than on the surface; at 6 m depth the 
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Table 3 
Soil type a v L z 
group m2/h mm/h m 
I 0, 0022 1,78 16,5 6,2 
ll 0, 0032 2,14 18,8 71,6 
lll 0, 0047 2,60 22,7 9,1 


coldest time is the beginning of July and the highest temperature occurs in the 
beginning of December, that is, approximately at times when the opposite extreme 
occurs on the soil surface, 

The temperature graphs in Figs. 4, 5 and 6 are based on the assumption that 
the soil surface is kept free of snow throughout the year, They approximate actual 
conditions also in winters with little snow, as can be seen from Figs, 7 and 8, re- 
presenting measured soil temperatures in till and clay soil in South-Pinland in 1957. 
Fig. 9 shows soil temperatures measured in the winter of 1956/57 in till and clay 
soil in South-Finland (10 km west of Helsinki), fn terms of depth from the soil 
surface. 

The families of curves in Figs. 4, 5 and 6 have been derived for homogeneous 
soil types; they are not quite accurately valid for stratified soils if the different 
layers have markedly different temperature conductivities. 

According to the investigations of Ad, Schmidt [19], the annual variation of 
temperature in homogeneous soil can be represented theoretically, with high accu- 
racy, when mean values of several years are concerned, In Canada (Saskatoon and 
Ottawa), D.C, Pearce [20] has carried out measurements of the annual variation 
of the temperatures in the ground; he found that already the temperatures calcu- 
lated by the first harmonic wave are well consistent with the measured values at 
2,4 m depth, the greatest deviation being 0,5°C, At 60 cm depth, already, the 
correspondence is almost equally good, i.e, about 1°C maximum deviation at 
all other times except in the spring, in April, when the thawing of the snow de- 
lays the heating of the ground, At this time the soil temperatures fall up to 4°C 
below the theoretical value. 

The depth of penetration of ground frost is dependent on the temperature cone 
ductivity of the soil and on the freezing index [% days). 

Fig. 10 shows the theoretical maximum variation of the frozen-ground limit depth 
under bare ground, the freezing of the soil moisture having been neglected, for 
three values of a in parallel with the freezing index, which has been computed as 
a mean for the above-mentioned ten-year period, 

The variation of the frozen ground limit in snow-covered ground, measured on 
a small moraine hill and on a field with clayey base soil in the previously men- 
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moraine hill in 1957. 


Fig- 7 Measured annual variation of temperatures at different depths in the snow-covered ground of a 
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Fig. 8 Measured annual variation of temperatures at different depths in a snow-covered clayey field 
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Fig. 9 Soil temperatures measured in till and clay soil in 1957, in depend- 
ence of depth 


tioned location (10 km west of Helsinki) in the winter of 1956/57 has also been 
plotted in Fig. 10, as well as the increase of the freezing index in the course of 
the winter. The snow cover of 30 cm average thickness has reduced the depth to 
which the ground froze to about one third of the theoretical value applying to bare 
ground, 
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Fig. 10 Theoretical limit depths of frozen ground (on the basis of Figs. 4, 
5 and 6); in the winter of 1956/57 for a = 0,0022 m#/h (Curve No, 1), 
a = 0,0032 m2/h (Curve No. 2) and a = 0,0047 m2/h (Curve No. 3). 


Measured limit depths of frozen ground in wll soil (Curve No, 4) and clay 

soil (Curve No. 5). Mean freezing index of the ten-year period 1945 « 1954 

(Curve No, 6), freezing index of 1956/57 (Curve No. 7) and thickness of 
snow cover in 1956/57 (Curve No, 8). 


Measurements of the depth to which the ground freezes have been carried out 
in Ottawa (Canada) in the winter of 1950/51 in various soil types in snow-covered 
and bare ground (Fig. i1) [2]. The figure also shows the variation in thickness 
of the snow cover and the freezing index. 
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Fig, 11 Variations of limit depth of frozen ground, thickness of snow cover 


and freezing index in Ottawa (Canada) in the winter of 1950/51. 
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Comparison of the frozen ground depth curve in Fig. 11 with the theoretical 
curve in Fig, 10 reveals that the frozen ground depth under bare ground is consid- 
erably less than the theoretical prediction, despite the higher freezing index. This 
is due to the fact that freezing of the water in the upper layers of the ground de- 
lays freezing in the autumn and thawing of the ice delays the heating of the soil 
in the spring, 

Equations have been derived for the mathematical calculation of the depth to 
which the ground freezes [15]. However, all such formulas are approximate be- 
cause the calculations are rendered inaccurate by the latent heat of the water and 
by the movements of water in the ground, More reliable results are obtained by 
estimating the frozen ground depth on the basis of the diagrams in Fig. 10 and 
Fig. 13 and of the results of measurement, 

Ground water produces changes in the temperature distribution in the ground if 
the ground water table lies above the depth to which the annual temperature var- 
iation penetrates. Water causes considerable changes in the temperature conduc- 
tivity of the soil. Moreover, flowing ground water acts as a heat-transporting me- 
dium. 

A shallow soil layer above a high-lying ground water table attains stationary 
conditions in a comparatively short time; for instance, this occurs within one 
month in a soil layer of 130 cm thickness if a = 0,0022 m2/h (Table 3). Ground 
water, as well as rock, has an equalizing influence upon the temperatures of the 
overlying soil layers. 
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4. THE QUANTITY OF HEAT STORED IN THE GROUND 


Knowledge of the heat quantities stored in the ground during the warm season is 
necessary, as they determine the temperature conditions in the ground during the 
winter, 

The quantity of heat stored in the ground at any time can be calculated if the 
annual temperature variation at different depths is known. 


Let us denote 
t=t, the annual mean temperature in the ground, 


t = t(z), the temperature at the depth z in the ground, 


The heat quantity stored in the ground is 


Figs. 12, 13 and 14 show the soil temperatures, in terms of depth, determined 
on the basis of Figs. 5, 6 and 7 on the Ist day of each month, By determining 
the areas of the curves with the planimeter, the corresponding values of the inte- 
gral (15) are obtained; multiplying them further with C [Mcal/m® °c] we find the 
values of Q. The annual variation of the heat stores in the ground are shown in 
Fig. 15 for three different values of the temperature conductivity. The heat quan- 
tities have been calculated with reference to the mean temperature of the previ- 
ously mentioned ten-year period, including depths up to 10 m. The maximum 
variation of the heat stored in the ground, or the total heat exchange of the ground 
has been shown in Table 4 for bare ground and for the different soil types. 


i 


Fig. 12 Temperatures at different depths in bare ground on the Ist of each 
month according to Fig. 4; a = 0,0022 m2/h. 


Table 4 
A Qmax 2min 
m2/h kcal/mh °C Mcal/m® °c Mcal/m? 
0, 0022 1,43 0, 650 27,17 
0, 0032 1,90 0, 594 30,0 
0,0047 2, 50 0, 532 34,1 


The function representing the annual variation of the soil surface temperature 
(formula (14) ) can also be used to calculate approximately, on the basis of the 
temperature gradient, the quantity of heat that has been absorbed in the soil 
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Fig. 13 Temperatures at different depths in bare ground on the ist of each 
month according to Fig, 5; a = 0,0032 m2/, 


t 
2 
T 
dt (16) 


z=0 


The heating of bare ground starts about the middle of March, and the heat 
stored in the ground is at its maximum in the beginning of September. When 


aT 
(34> = is calculated from equation (13), taking into account only its first 


cosine term, which is graphically represented by Curve No. 3 in Fig, 3, and when 
this is substituted in equation (16), one finds by integration 
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Fig. 14 Temperatures at different depths in bare ground on the Ist of each 
month according to Fig. 6; a = 0,0047 m2/h, 


Substituting in (17), 


T = 12,2°C tt, = 15 days 
to = 365days tg = 250 days 


the heat quantity stored in the ground is obtained for the different soil types (Ta- 
ble 5). The heat quantities are 10 -- 15 % lower than the values in Table 4, which 
were calculated by a more accurate procedure. 

According to different investigators, the total heat exchange of the ground is 
1530 Mcal/m?, depending on conditions, calculated as the difference of the 
maximum and minimum heat quantities contained in the ground, when the heat 
capacity varies between 0,3 and 0,5 Mcal/m? °C, Most nearly equivalent to con- 
ditions in South-Finland are the measurements performed by Luboslavsky [ 22] in 
the vicinity of Leningrad. His calculations are based on mean temperature values 
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a A Q 
m2/h kcal/mh °C Mcal/m2 
0,0022 25, 1 
0, 0032 1,90 21,6 
0, 0047 2, 50 30,0 


measured during 15 years at depths between 0 and 160 cm. According to Lubo- 
slavsky, sandy soil, with the specific heat 0,57 Mcal/m® °C, has a total heat 
exchange of 


Qmax ~ Qmin = 26,51 Mcal/m? (18) 


The greater part of the heat is stored in the upper layers of the ground, a con- 
siderable proportion being bound as latent heat in the ground layer which freezes 
in the winter. In till soil, for instance, 73 % of the total exchange quantity are 
stored in the 0+ 160 cm layer according to the previous theoretical calculation. 
Luboslavsky reports for sandy soil the corresponding figure of 61 % The soil type, 
especially its temperature conductivity, has an essential influence upon the heat 
quantity stored by the soil. High ground water table increases the heatestoring ca- 
pacity of the ground. 

Plant cover in the summer hinders the heating of the ground and snow cover in 
the winter hinders its cooling. According to Luboslavsky, the difference of the 
maximum and minimum stored heat quantities in a soil layer of 160 cm thickness 
is about 33 % less for ground covered by vegetation and by snow than for bare 
ground, Fig. 16 shows the heat quantities stored in bare and covered ground, by 
the seasons, according to Luboslavsky’s results. Between September Ist and February 
28th, ground protected by snow cover gives off only 46 % of the heat quantity es- 
caping during this time from bare ground. In the summer, again, snow cover de- 
lays the heating of the ground, as it reflects part of the incident thermal energy 
and part of it is consumed for the thawing of snow; consequently, the temperature 
in snow-covered ground begins to rise 14 month later than that of bare ground. 
Upon thawing of the snow, the protection provided by vegetation is insignificant 
yet, and the ground is therefore warmed up rapidly. During the actual heating 
period about 79 % of the heat quantity correspondingly taken up by bare ground is 
stored in covered ground, In the autumn, when the protective effect of vegetation 
has become less, covered and bare ground cool off at about the same rate. 

According to the theoretical calculation (Fig. 15), the maximum heat quantity 
is stored in the ground at the beginning of September when the heating period of 
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houses begins. The heat stores of the ground are at minimum about the 20th of 
March. During the interval the ground has given off the entire stored heat quan- 
tity (Table 4), which benefits the buildings. 

The heat flow from the ground to the surface, for bare soil, is obtained from 
Fig. 15 for the different soil type groups as the derivative with respect to time of 
the function representing the stored heat quantity at any given time. The heat 
flow is highest, 8 «11 kcal/m? h, at the turn of the year. In snow-covered ground, 
the corresponding values are obtained from the temperature gradient of the meas- 
ured field at the surface. For instance, ot = 3 °C/m (Fig. 9) and when the ther- 
mal conductivity isA = 1,9 kcal/mh °C, the upward flow of heat is on an aver- 
age 5,7 kcal/m* h for snow=cover varying between 40 and 60 cm. According to 
measurements carried out by L.W. Cold [23] in Ottawa, Canada (@ = 45°24’ 
N.la., Y = 15 943°), the total heat quantity of 24500 kcal/m?2 flowed out from 
the ground during the time November 17th, 1956 to March 3rd, 1957 at a snow 
depth of 20 cm on an average. This is equivalent to an average heat flow of 
9 kcal/m? h, 
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